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1^ ■ Abstract 

. We compare results of 2 — > 1 /c^-factorization approach with Kwiecihski unintegrated parton 

o . 

distributions and the standard cohinear factorization approach at RHIC and slightly smaller en- 
' ergies. Our approach contains only one free parameter responsible for internal parton motion in 

nucleons. In contrast to recent works in the literature our fct-factorization approach includes also 
CJ \ quark degrees of freedom in addition to purely gluonic terms. Both mid and forward rapidity 

^ ' regions are considered. We discuss uncertainties due to fragmentation functions. In general, the 

/ct-factorization approach gives a better description of the ~ 1 - 4 GeV region both at mid and 
forward rapidity regions. Our approach leads to asymmetry in the production of vr"*" and vr", very 
similar to the one observed very recently by the BRAHMS collaboration. 

PACS numbers: 12.38.Bx, 13.85.Hd, 13.85.Ni 
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I. INTRODUCTION 

In order to extend the applicability of the collinear-factorization approach (see e.g. 0, 
to jet and/or meson production to small transverse momenta it was proposed 
io add an extra Gaussian distribution in the transverse momentum of colliding partons 
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12| |. In addition to assuming the parameter of the Gaussian distribution 
to be X- independent, one has to substitute t —>■ t — fi^, u ^ u — fi"^ and s —>■ s + 2yU^ in 
order to avoid infinities in the matrix element squared contributing in this approach at finite 
transverse momenta. This procedure is an artifact of the on-shell approximation (on-shell, 
non-coUinear partons) used. The result of such a procedure depends on the value of the 
parameter /x^, especially at low transverse momenta. A part of the arbitrariness can be 
avoided in the fcf-factorization approach (see e.g. |13'] ) . 

The fc^-factorization approach was used recently to study both the rapidity and trans- 



14 llq with only 



verse momentum distributions of particles produced at RHIC energies 
gluon degrees of freedom taken into account. Recently we have supplemented the mecha- 
nisms inclu ding the gluon degrees of freedom with mechanisms including the quark degrees 



of freedom 



13|. The latter are very important at very forward [r] ^ 0) and very backward 



{f] ^ 0) rapidity regions. Our approach makes use of the recently developed Kwieciiiski 



unintegrated parton distributions 
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18l |. In contrast to other approaches in the liter- 



ature the formalism takes into account the x-dependent radiative effect of kt - broadening 
of "initial" parton distributions. The formalism of the Kwiecihski UPDF is adequate in the 
region of not too small longitudinal momentum fractions i.e. at not too high energies. We 
lave applied these unintegrated parton distributions to gauge boson IiqI.^ standard Higgs 
2ol | hadroproduction as well as for charm-anticharm photoproduction [21]. Very rece ntly 
we have applied this formalism also to the description of the SPS pion production data jl3| . 
In the present analysis we shall apply it to the pion hadroproduction at somewhat larger 
energies, up to the RHIC energy W = 200 GeV. 

II. INCLUSIVE CROSS SECTIONS FOR PARTONS 

The approach proposed by Kwiecihski is very convenient to introduce the nonperturba- 
tive effects like internal (nonperturbative) transverse momentum distributions of partons in 
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nucleons. It seems reasonable, at least in the first approximation, to include the nonpertur- 
bative effects in the factorizable way 



Mx,b,fi') = fr\x,b,f,')-Fr{b) 



(1) 



where the index i denotes either gluons or quarks or antiquarks. The form factor responsible 
for the nonperturbative effects must be normalized such that 



= 0) = 1 . 



(2) 



In the following, for simplicity, we use a flavour and x-independent form factor 



F^P{b) = F^P{b) = F^P{b) = = exp - 



Abl 



(3) 



which describes the nonperturbative effects. The Gaussian form factor in b means also a 
Gaussian initial momentum distribution oc exp(— A;^6q) (Fourier transform of a Gaussian 
function is a Gaussian function). Gaussian form factor is often used to correct coUinear 
pQCD calculations for the so-called internal momenta. Other functional forms in b are also 
possible. 

In the fcj-factorization approach usually the gg ^ g fusion mechanism is included only 
Q. In Ref.Q we have included two other leading-order diagrams which involve quark 
degrees of freedom. They are important in the so-called fragmentation region 

Q- The 

momentum-space formulae for all the processes included read: 
for diagram A (gg^g): 



da' 



l6Nc 1 



dydPpt ^1 — Ipl 

j Ols{^^) fg/liXl,Kl,fX^) fg/2{X2, (1^) 
+ K2- Pt) d^Kid^K2 , 



(4) 



for diagram Bi (q/ g q/): 



-)- 

dyd^Pt iV2 - 1 \<d) pj 



E 



) fqj/li.Xl,lii,ll ) /<;/2(a;2,K2,/X ) 



+ K2 - Pt) d^Kid^K2 



(5) 
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for diagram B2 (g q/ — ^ q/): 

da^^ _ 16N^ /4\ 1 



/ 

+ K,2— Pt) Cpl^l(f'K2 ■ (6) 



These seemingly 4-dimensiona 
suitable change of variables |l5| 



integrals can be written as 2-dimensional integrals after a 



j ... 5^'\K^ + K2-Pt) d'K,d'K2 = I ... ^ . (7) 

The integrands of these "reduced" 2-dimensional integrals in qt = ki — K2 are generally 
smooth functions of qt and corresponding azimuthal angle (pq^. In Eqs.(jl)), © and (jHI) the 
longitudinal momentum fractions 

= ^"^ t- exp(±g) . (8) 

n 

where rrix is the effective mass of the parton. This is important only at — > [13|. 

The sums in (jSJ and © run over both quarks and antiquarks. The argument of the 
running coupling constant above was not specified explicitly yet. In principle, it can be 

or a combination of pf, and /t^. In the standard transverse momentum representation 
it is reasonable to assume fl^ = min(p^, nf, K2) (see e.g. In the region of very small pt 

usually p^ < k\, and ^2 = is a good approximation. 

Assuming for simplicity that Vt^ = Q'^{p'f) or pf (function of transverse momentum squared 
of the "produced" parton, or simply transverse momentum squared) and taking the following 
representation of the 5 function 

S^'^\ki + K2- Pt) = j d% exp {ki + K2 - Pt)b , (9) 

the formulae Q, (0) and © can be written in the equivalent way in terms of parton 
distributions in the space conjugated to the transverse momentum. The corresponding 
formulae read: 



for diagram A: 



da^ IQN^ 1 , ON 

T<^s{Pt 



dyd^Pt N^-lpr 
j fg/i{xi,b,fx'^)fg/2ix2,b,n^)Jo{ptb)27rbdb, (10) 



for diagram Bi 



dyd^Pt ~ - 1 V9 J Pt 
^ J fgf/iixi,b,(x'^)fg/2ix2,b,fx'^)Jo{ptb)2nbdb, (11) 



for diagram B2: 



77 -o^siPt] 



da^^ _ IQN, _ _ 
dyd^Pt ~ - 1 [9) ft 

^ / fg/i{xi,b,n^) fg^./2{x2,b,^^)Jo{ptb) 2TTbdb . (12) 
/ 

These are 1-dimensional integrals. The technical price one has to pay is that now the 
integrands are strongly oscillating functions of the impact factor, especially for large pt- 
The formulae (jlUp . and (|12p are very convenient to directly use the solutions of the 
Kwiecihski equations discussed in the previous section. 

When exten ding running ag to the region of small scales we use a parameter-free analytic 



model from ref. 



221. 



III. FROM PARTONS TO HADRONS 



In Ref . |14| it was assumed, based on the concept of local parton-hadron duality, that 
the rapidity distribution of particles is identical to the rapidity distribution of gluons. In 
the present approach we follow a different approach which makes use of phenomenological 
fragmentation functions (FF's). In the following we assume 9^ = 9g. This is equivalent to 
Vh = Vg = Vgi where rj^ and rjg are hadron and gluon pseudorapitity, respectively. Then 

Ug = arsinh ( sinh yh ] , (13) 

V Pt,h I 



where the transverse mass mt^h = y ^'h, + pIk- order to introduce phenomenological FF's 
one has to define a new kinematical variable. In accord with e^e" and ep collisions we define 
a quantity z by the equation = zEg. This leads to the relation 



Pt,g = —J{mt,h,yh) , (14) 



z 



where the jacobian Jimt^hjUh) reads 



-1/2 



J(m,,, y,) = 1 - — ^ . (15) 

m^^cosh yh 



Now we can write a given-type parton contribution to the single particle distribution in 
terms of a parton (gluon, quark, antiquark) distribution as follows 

daP{r]h,Pt,h) 



dcTiyp,Pt,p) 



zpt,p 
J 



dypd^Pt-p 



(16) 



Please note that this is not an invariant cross section. The invariant cross section can be 
obtained via suitable variable transformation 

dcr^{yh,Pt,h) _ ( d{yh,Pt,h)\~^ daP{yh,Pt,h) 



where 



dyhd'^Pt,h 



yh = 2 



d{Vh,Pt,h] 



drihd'^Pt,h 



(17) 



p'ih 



+ sinh rjh + sinh rjh 



nii +p, . o 

' + smh ijh — smh ijh 



p'i,h 



(18) 



Making use of the 6 function in (fT^ the inclusive distributions of hadrons (pions, kaons, 
etc.) are obtained through a convolution of inclusive distributions of partons and flavour- 
dependent fragmentation functions 

da{r]h,pt,h) 



drjhd^pu. 



dz— 
z^ 



/=-3 
3 

/=-3 



dn^^ 

I 2\ 'if9-'gf^y'if'^t'if' 



dygd^Pt^g 

XVqf^PtA 



dyqjd^Pt^q 
d(Tgg^^.^q^.{yqj,Pt,qf) 



dygjd'^Pt,, 



yq=-nh 



(19) 



One dimensional distributions of hadrons can be obtained through the integration over the 
other variable. For example the pseudorapidity distribution is 



daijh 
drjh 



d^Pt,h 



da{r]h,pt,h) 
drjhd^Pt^h 



(20) 



There are a few sets of fragmentation functions available in the literature (see e.g. 
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IV. RESULTS 



23 



24 



In the present application different fragmentation functions from the hterature 
will be used. All of them were obtained in global fits to the e^e~ data. In the present paper 
we shall show related uncertainties for pion hadropro duct ion. 

Before we go to the RHIC data we wish to look at some lower-energy data measured at 
ISR 27]. It was pointed out recently [5J that the standard collinear approach is not able 
to describe very forward production of vr''. In Fig^we present invariant cross section as a 
function of the Feynman for W = 23.3 GeV at two different laboratory angles 9 = 15°, 
22°. The collinear factorization result (left panel) lies well below the data independent of 
what fragmentation functions are used. The results of the fc^-factorization approach describe 
the experimental data much better. In Fig|21we present similar results for somewhat larger 
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FIG. 1: Invariant cross section as a function of at W = 23.3 GeV for different angles. The 
standard collinear result is shown in panel (a) and result of our approach in panel (b). The 
experimental data are from [2^ 

energy W = 52.8 GeV and 9 = 5°, 10°. The situation here is very similar to that at the 
lower energy. 

Clearly, the standard collinear approach fails badly in the very forward region. The kt- 
factorization is much better but still some contribution is missing at large Feynman xp and 
at small transverse momenta. However, one should remember that other processes such as 
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FIG. 2: Invariant cross section as a function of xp at W = 52.8 GeV for different angles. The 
standard collinear result is shown in panel (a) and result of our approach in panel (b). The 
experimental data are from [2^ . 



pion stripping j36| and/or diffractive production of nucleon resonances and their subsequent 
decay may play an important role here. This requires a separate analysis which goes beyond 
the scope of the present paper. 

Let us return to the midrapidity region. The PHENIX collaboration has measured in- 
variant cross section as a function of the vr^ transverse momentum at W = 200 GeV in a 
very narrow interval of pseudorapidity rj = 0.0 ± 0.15. 

In FigBlwe show our full result (diagrams A, Bi and B2 Q|) for different fragmentation 
functions 2^ . 2^. In this calculation ho = 0.5 GeV^^ was used. This is the optimal value 
of the parameter for gauge boson production jl^. In Fig 01 we show the dependence on the 
value of the parameter Bq. Having in view the uncertainties in the fragmentation functions 
and in the parameter bo, which is responsible for nonperturbative effects like e.g. parton 
Fermi motion, our fc^-factorization result describes the data very well and in a quite broad 
range of transverse momentum. In FigOlwe show individual contributions of gluon and quark 
components. In contrast to the standard beliefs (see e.g. the quark contributions are 
only slightly smaller then the gluon ones and certainly not negligible, especially at larger 
transverse momenta. 
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FIG. 3: Invariant cross section for vr'^ production as a function of pion transverse momentum 

at W = 200 GeV and ri = 0.0. The fc-factorization results are shown for different values of the 

□ 

parameter bo. The experimental data of the PHENIX collaboration are from |28l |. 

Let us come now to charged particle distributions. In FiglHlwe show the results of the 
collinear approach for different values of hadron pseudorapidity for the Kretzer fragmentation 
functions 2J]. The results of the calculations are compared to the BRAHMS collaboration 
data 29]. The theoretical calculations are for the charged pions. The data are not exactly 
for pions but include all charged hadrons, e.g. protons as well. There is a large factor of 
disagreement between the calculations and the data. This factor seems to increase when 
going from mid to forward rapidities which is not compatible with the NLO collinear factor- 
ization approach where the so-called K-factor is almost independent of rapidity. In Fig|7|we 
)resent result of our calculations for Bq = 0.5 GeV and the Kretzer fragmentation functions 
2^ . While at large (pseudo)rapidities {rj = 2.2, 3.2) our results (negative pions) are quite 
compatible with the experimental result for negative hadrons, there seems to be a missing 
contribution at more central (pseudo)rapidities [r] = 0.0, 1.0), i.e. in the case when the sum 
of positively and negatively charged hadrons is measured. At present it is not clear to us if 
the missing strength is due to protons and/or positively charged kaons. In Fig|Hl individual 
contributions from our fct-factorization approach are shown. They correspond to diagrams 
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FIG. 4: Invariant cross section for vr^ production as a function of pion transverse momentum at W 
= 200 GeV and r] = 0.0. The A:j-factorization results are shown for different sets of fragmentation 



functions. The experimental data of the PHENIX collaboration are from 



2S|. 



A, Bi and B2 in Fig.l of Ref . [isj] . While at = the gg ^ g contribution dominates up 
to pt ~ 3.5 GeV at rj = 2.2 the Bi contribution is larger than the gg ^ g one already at 
Pt ~ 1.5 GeV and at 77 = 3.2 above ~ 1 GeV. This is the Bi contribution which provides 
a good description of the BRAHMS data at forward {1] ^ 2) rapidities. 

In Figiniwe compare results of our calculations with the old proton-antiproton charged 
hadron data 31] for rj = 0. The dependence on fragmentation functions is shown in the left 
panel. In the right panel we show the dependence on the parameter Bq for the Kretzer frag- 
mentation functions. The figure looks very similar to the figure with the BRAHMS collabora- 
tion data for 1] = 0, which simply reflects consistency of the proton-proton (BRAHMS) data 
and proton-antiproton (UAl) data. Within the approximations used in our approach the 
charged pion inclusive cross section is identical for the proton-proton and proton-antiproton 
collisions, provided the energy is the same, which is the case for the BRAHMS and UAl 
colaboration data. 

Let us concentrate now at the very forward region of the phase space. Recently the STAR 
collaboration has published large-rapidity, intermediate-transverse-momentum data for 
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FIG. 5: Invariant cross section for vr production as a function of pion transverse momentum at W 
= 200 GeV and ij = 0.0. Individual contributions of gluon and quark fragmentation are shown for 
different fragmentation functions. The experimental data of the PHENIX collaboration are from 
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FIG. 6: Invariant cross section for charged particle production for different values of particle 
pseudorapidity at W = 200 GeV. The lines represent results obtained within standard collinear 

n 

factorization approach and Kretzer fragmentation functions 24]. The BRAHMS collaboration 



experimental data [29| are shown by the solid circles. 



7r° production. In Fig^J we compare the results of the colhnear (left panel) and our kt- 
factorization approach (right panel) calculated with the Kretzer fragmentation functions. 
While the collinear approach underestimates the STAR experimental data by a factor of 
about 3, the fcrfactorization approach is almost consistent with the data, especially at 
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FIG. 7: Invariant cross section for charged particle production for different values of particle 
pseudorapidity at W = 200 GeV. The lines represent results obtained within our fct-factorization 
approach for different values of parameter bo = 0.25, 0.5, 1.0 GeV~^ and Kretzer fragmentation 
functions 24]. The BRAHMS collaboration experimental data j^] are shown by the solid circles. 



smaller pion energies, i.e. at not too high transverse momenta. The situation would change 
somewhat with different set of fragmentation functions. 

In our approach the diagrams with quark degrees of freedom {Bi and B2) lead to an 
asymmetry in 7r+ and vr" production. Recent results of the BRAHMS collaboration for the 
transverse momentum integrated cross section da/dy show the — n~ asymmetry already 

i i 

above y=2 [32]. As an example in Fig^2 we show the ratio of the corresponding cross 
sections for vr" and 7r+ production for the STAR kinematics. A huge deviation from the 
unity can be observed. The /cj-factorization approach ratio (thick solid line) is somewhat 
smaller than the corresponding ratio in the coUinear approach (dashed line). 

During the preparation of this manuscript the BRAHMS collaboration has presented 
the first preliminary data for identified vr"*" and vr" jsSj. In Fig^J we show transverse 
momentum dependence of the ratio 71^/71^ for rapidity y = 3.1. Only statistical error bars 
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FIG. 8: Invariant cross section for charged particle production for different values of particle 
pseudorapidity at W = 200 GeV. The lines show individual contributions within our /ct-factorization 



approach. Here Kretzer fragmentation functions 
experimental data 3] are shown by the solid circles. 



2J] were used. The BRAHMS collaboration 



are shown. Although the experimental ratio is not completely monotonous it clearly shows 
a deviation from unity which is an experimental evidence that gluon degrees of freedom are 
not sufficient to describe the production of hadrons. We describe the main trend of deviation 
of the ratio from unity relatively well except in the region of very small values of transverse 
momenta, js?! Of course such a ratio depends on the details of the fragmentation functions 
and in particular on their flavour decomposition which as discussed above is very difficult 
to obtain from the e^e~ scattering alone. We hope that in the near future the BRAHMS 
collaboration will be able to scan the ti^/ti^ ratio as a function of (pseudo)rapidity and 
transverse momentum in order to identify the contributions with quark degrees of freedom 
and perhaps to extract flavour-dependent fragmentation functions. 

In nuclear collisions the 7r+ — tt" asymmetry is weakened by the presence of the pp, nn, 
pn and np subcollisions. Due to isospin symmetry relation, this leads to equal yield of 7r+ 
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FIG. 9: Invariant cross section for charged particle production for different values of particle 
pseudorapidity at W = 200 GeV. The uncertainties on fragmentation function are discussed in 
the left panel and the dependence on the 60 parameter in right panel. The UAl collaboration 



experimental data are shown by the solid circles. 

and vr" for collisions of isospin symmetric nuclei. For collisions of heavy nuclei there is a 
small nonzero effect due to the excess of neutrons over protons. The asymmetry can be, 



however, quite sizeable in peripheral collisions [3^ due to neutron skin effects. The sign of 
the nuclear asymmetry is then reversed as compared to the proton-proton collisions. The 
small asymmetry in nuclear collisions was recently mistakenly interpreted as a dominance of 
purely gluonic effects even at large rapidities. Our calculation actually shows that the quark 
terms dominate at very forward/backward rapidity regions. Although our calculation is for 
elementary proton-proton collisions only it puts into question some recent nuclear color glass 
condensate calculations based on gluon degrees of freedom only. 

V. CONCLUSIONS 

We have shown that the formalism recently developed by us and based on unintegrated 
parton distributions which fulfill the so-called Kwiecihski evolution equations provides a rea- 
sonable description of the recent experimental data of the PHENIX, BRAHMS and STAR 
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FIG. 10: Invariant cross section for neutral pion production as a function of pion energy at W = 
200 GeV for rj = 3.8 for collinear (left panel) and our fet-factorization (right panel) approach. The 
STAR collaboration data are taken from Ref.j^^. 



collaborations at RHIC. The description is particularly good in the region of intermediate 
transverse momenta of pions pt^h ~ 1 - 4 GeV. In comparison, the standard collinear fac- 
torization approach gives results which are by a factor of 3 - 7 lower than the experimental 
data, depending on particle rapidity. We have found a rather strong dependence on the set 
of fragmentation functions used in the calculation. 

Inclusion of diagrams with quark degrees of freedom leads to 7r+ — tt~ asymmetry. The 
preliminary BRAHMS data provide evidence for such an asymmetry. A dedicated measure- 
ment of the 7r+ — TT^ asymmetry in forward and backward region as a function of transverse 
momentum would be a good test of the present approach and perhaps could be used to con- 
strain better the gluon-to-pion and quark-to-pion fragmentation functions which extraction 
in e~^e~ collisions is ambigous to a large extent. 
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FIG. 11: The ratio of the vr to vr"'" cross sections as a function of pion energy for r/ = 3.8 for 
cohinear (dashed) and our /cj-factorization (thick sohd) approach. 
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